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ABSTRACT 

We investigate the evolution of angular momentum in simulations of galaxy forma- 
tion in a cold dark matter universe. We analyse two model galaxies generated in the N- 
body/hydrodynamic simulations of Okamoto et al. Starting from identical initial conditions, 
but using different assumptions for the baryonic physics, one of the simulations produced a 
bulge-dominated galaxy and the other one a disc-dominated galaxy. The main difference is the 
treatment of star formation and feedback, both of which were designed to be more efficient in 
the disc-dominated object. We find that the specific angular momentum of the disc-dominated 
galaxy tracks the evolution of the angular momentum of the dark matter halo very closely: the 
angular momentum grows as predicted by linear theory until the epoch of maximum expan- 
sion and remains constant thereafter. By contrast, the evolution of the angular momentum of 
the bulge-dominated galaxy resembles that of the central, most bound halo material: it also 
grows at first according to linear theory, but 90% of it is rapidly lost as pre-galactic fragments, 
into which gas had cooled efficiently, merge, transferring their orbital angular momentum to 
the outer halo by tidal effects. The disc-dominated galaxy avoids this fate because the strong 
feedback reheats the gas which accumulates in an extended hot reservoir and only begins to 
cool once the merging activity has subsided. Our analysis lends strong support to the classical 
theory of disc formation whereby tidally torqued gas is accreted into the centre of the halo 
conserving its angular momentum. 
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1 INTRODUCTION 

In the standard theory of galaxy formation in the cold dark mat- 
ter (CDM) cosmogony, disc galaxies form when rotating gas cools 
inside a dark matter halo and fragments into stars. Analytic and 
semi-analytic calculations of this process start by assuming that 
tidal torques at early times impart the same specific angular mo- 
mentum to the gas and the halo and that the angular momentum of 
the gas is conserved during the collapse (Fall & Efstathiou 1980l; 
IMo et alJll998l : ICole et al. ll200d : lFirmani & Avila-Reeseii200ClO . It 
is unclear whether this assumption holds in gasdynamic simula- 
tions of galaxy formation. 

The first attempts to simulate the formation of a spiral galaxy 
from cold dark matter initial conditions generally failed, produc- 
ing objects with overly centrally concentrated distributions of gas 
and stars. It became immediately apparent that the root cause of 
this problem was a net transfer of angular momentum from the 
baryons to the dark matter halo during the aggregation of the galaxy 
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through mergers jNavarro & Benj|l99ll : iNavarro & White! 1 19941: 
iNavarro. Frenk. & Whitell995h . This is known as the "angular mo- 
mentum problem". It was suspected from the start that its solution 
was likely to involve feedback processes that would regulate the 
supply of gas to the galaxy. 

More recent simulations within the CDM framework have 
produ ced more promising disc galaxies. iThacker & CouchmanI 
1 120011) obtained a reasonably realistic disc at 2 = 0.5 (when 
their simulation stopped) by ass uming that gas cooling is stro ngly 
suppressed by feedback effects. ISteinmetz & Navarrq ( l2002h . us- 
ing different prescriptions for star formation and feedback, found 
that a broad ran ge of galax y morp hologies could be produced. 
In related work, lAbadi et al.f ( 120031) obtained a disc galaxy with 
a realistic density profile but no t enou gh angular momentum. 
ISommer-Larsen. Gotz. & Portinaril ( |2003|) found that considerably 
larger energy feedback rates at early times than previously em- 
ployed were required to ameliorate the angular momentum prob- 
lem. They achieved this by distinguishing between early and late 
star formation modes, assuming a higher star formation efficiency 
and stronger feedback in the early mode. They were able to gen- 
erate a variety of morphological types, including discs whose an- 



2 J. Zavala, T. Okamoto, and Carlos S. Frenk 



gular momentum, however, was still a factor of two smaller than 
obse rved for real disc galaxi es. 

iGovemato et alj | |200 4') produced a disc galaxy without re- 
quiring strong feedback and ascribed the angular momentum prob- 
lem to a lack of num erical resolution. However, subsequently, 
iGovemato et al.l ^20071) were forced to invoke strong feedback in 
order for their simulated discs to lie closer to the Tully-Fisher re- 
lation. While resolution is c learly an important and complicated 
issue, lOkamoto et al. I ( I2OO3I) argued that if the hot and cold gas 
components are decoupled from each other, then the effects of 
resolution in th e curre nt generation of simulations are reduced. 
[Robertson et al.l ( l2004t) adopted a multiphase model for the star- 
forming interstellar medium (ISM) which stabilises gaseous discs 
against the Toomre instability, and produced a galaxy with an expo- 
nential surface brightness profile but insufficient angular momen- 
tum in the inner parts. They concluded that the fragmentation of 
gaseous discs caused by the Toomre instability is the main cause of 
the angular momentum loss in simulations that usually assume an 
isot hermal ISM at a tempe rature of ~ 10'* K. 

lOkamoto etalj ( l2005b also adopted a multiphase model for the 
ISM and assumed a top-heavy initial mass function (IMF) for stars 
formed in starbursts. T he latter assumptio n was motivated by the 
semi-analytic model o f iBaugh et all ( l2005h which requires a large 
relative abundance of massive stars in order to match the number 
counts of bright submillimeter galaxies. A top-heavy IMF is fur- 
ther supported by the good agreement of this model with the metal 
abundances measured in the intr acluster medium tNagashima et al.l 
l2005ah and in elliptical galaxies jNagashima et alj2005bh. By vary- 
ing th e criteria for starbursts in their simulations. lOkamoto et al.l 
( l2005h were able to produce galaxies with a variety of morphologi- 
cal types, from ellipticals to spirals, starting from exactly the same 
initial conditions. They showed that disc galaxies can form in halos 
that have a relatively quiet merger history if the early collapse of 
baryons is inhibited by strong feedback, in this case generated by 
the evolution of a stellar population with a top-heavy IMF. 

In this paper, w e analyse two of the galaxies simulated by 
lOkamoto etai] ( l2005h . one which ended up as a bulge-dominated 
galaxy and another that formed an extended stellar disc. We trace 
the time evolution of the angular momentum of the dark matter 
and the gas, and we distinguish between the evolution of the in- 
ner part of the dark halo and the halo as a whole. We show ex- 
plicitly that mergers transfer angular momentum from the inner 
parts of the system to the outer halo. This process is one aspect 
of the redistribution of angular momentum during virialization in- 
vestigated in N-body simul ations by D'Onghia et al. (2006) and 
iD'Onghia & Navarrol J2007h . Thus, the spin of the central part of 
the final system contains a fossil record of the merger history of the 
object. 

This paper is organised as follows. The simulations are briefly 
described in Section 2. The method we use to follow the angular 
momentum of the different components is explained in Section 3 
where we also present our results. We conclude in section 4. 



2 THE SIMULATIONS 

Details of the s imula tions that we analyse here may be found in 
lOkamoto et all ( l2005h . The initial conditions were extracted from 
a cosmological N-body simulation (of a cubical region of side 
35.325ft~^ Mpc) and correspond to a dark matter halo of present 
day virial mass Mvir — 1.2 x 10^^ Mq which has a quiet merg- 
ing history since z ^ 1. The simulations were carried out using 





— 







-20 -10 10 20 - 20 -10 10 20 

« (»pcA> « C»P<:/I') 



Figure 1. Edge-on (top) and face-on (bottom) views of the baryonic com- 
ponent of the bulge-dominated (left) and disc-dominated (right) galaxies. 
Stars and gas with pg ^ 7 X lO^^^grcm^'' within 10% of the virial 
radius make up the baryonic component. The brightness indicates the pro- 
jected mass density. 



the N-body/SPH code GADGET-2 l lSDringelll2005l) . They include a 
multiphase description of the star-forming gas and make use of the 
'phase-decoupling' technique introduced by Okamoto et al. ( 200 J 
to suppress the spurious angular momentum transfer from cold disc 
gas to ambient hot halo gas which would otherwise arise from the 
very nature of SPH. 

One of the two simulations that we analyse is the "no-burst" 
model which has a sta ndard prescriptio n for star formation, tuned 
to reproduce the local lKennicuttl ( ll998t) relation between star for- 
mation rate and gas density. This simulation yielded a galaxy with 
a large bulge and a small disc, to which we will henceforth refer as 
the "bulge-dominated" galaxy (left panels in Fig.[T|l. 

The second simulation is the "shock-burst" model in which 
the star formation efficiency is higher than average in gas that 
has recently been shock-heated. The stellar populations that form 
from this gas in a starburst are assumed to have a top-heavy IMF. 
Thus, in this model there is additional supernova heating avail- 
able when galaxies merge. The outcome of this simulation was a 
galaxy with a large disc and a small bulge - the "disc-dominated" 
galaxy (right panels in Fig. [TJ. The (kinematically identified) disc 
in this galaxy contains 48% of the stellar mass and 84% of the 
B-band light, making this a good candidate for a late-type spiral. 
However, the surface brightness distribution of the disc is consid- 
era bly more extended tha n is commonly observed for such galax- 
ies ('Okam oto et"ai]|2005l) . Interestingly, the satellites of the main 
galaxy in this model provide a good match to the properties of the 
brigh test Milky Way satell ites dubeskind et all 20071) . 

lOkamoto et alj ( |2005|) presented a third simulation in their pa- 
per, the "density-burst" model, which led also to a bulge-dominated 
object. We do not discuss our analysis of this object in this paper 
because this galaxy lost most of its baryons in strong galactic winds 
and is therefore not particularly suitable for studying the evolution 
of the angular momentum of the baryonic component. Neverthe- 
less, we have checked that the results for this object are qualita- 
tively similar to those of the "bulge-dominated" galaxy. 
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Figure 2. Upper panels: the fraction of the final dark halo mass that is contained in progenitor halos with particle number. A^dm ^ 10> at redshift z for the 
bulge-dominated (left) and the disc-dominated (right) galaxies. Lower panels: baryon fractions in the progenitor halos. The thick solid lines give the total mass 
fraction of baryons in progenitor halos with ^ 10. The dotted line shows the contribution to this fraction from cold baryons defined to be stars and 

dense gas with pg ^ 7 X 10^^'^gr cm^'*; the dashed line shows the confiibution of all other baryons to which we refer as warm/hot gas. The thin line marks 
the cosmic baryon fraction, Q^,/Qo. The left panel corresponds to the bulge-dominated galaxy and the right panel to the disc-dominated galaxy. 



In our analysis of angular momentum evolution, we distin- 
guish three (Lagrangian) components for each galaxy at 2: = 0. 
The centre of the system is taken to be the density peak of the stel- 
lar distribution. The first component is the dark matter halo which 
we define as the dark matter mass contained within the virial ra- 
dius, r„ . We take r„ to be the r adius within whic h the mean density 
is 100 times the critical value teke. Cole, AFren k 1996). For both 
galaxies, the halo mass is ~ 10^^h~^ Mq and the virial radius is 
^ 210 /i~^kpc. The second component we consider is the inner 
Lagrangian halo which we define as the 10% most bound dark mat- 
ter particles at redshift z = 0. The third component is the "galaxy" 
which we define as the stars and dense gas (p ^ 7x lO^^^gcm"^) 
contained within O.lri,. In the next section, we will follow the evo- 
lution of these three components. 

Fig.[T]shows edge-on and face-on views of the surface density 
of the two galaxies. The bulge-dominated galaxy (left) has a cen- 
trally concentrated spheroidal bulge, a small disc and an extended, 
nearly spherical stellar halo. The disc-dominated galaxy (right) has 
a small central bulge, an extended disc and a spheroidal stellar halo. 

In Fig.|2]we highlight relevant aspects of the global evolution 
of the Lagrangian region associated with the final dark matter halo. 
We identify the progenitors of the final halo by applying, at each 



redshift, the friends-of-friends grouping algorithm of lOavis et al.l 
( 1985) to the dark matter particles that end up in the halo at z = 0. 
We employ a linking length, lunk — b(z){l), where {1} is the mean 
dark matter interparticle separation and the linking parameter, b{z), 
is scaled from the canonical value of 0.2 for an Einstein-de-Sitter 
universe according to 

where Av (z) is the virial overdensity (relative to the mean) at each 
redshift. As expected, the evolution of the dark matter component 
is very similar for the two simulations. The halo builds up rapidly 
at first and, by 2 = 4, half of the final halo mass is already in place 
in progenitors of more than 10 particles each. Large mergers are 
manifest as rapid increases in mass, such as that visible at 2 ~ 1.4. 

The evolution of the baryons in the collapsed regions is very 
different in the two simulations. Initially, all the baryons are as- 
sumed to be cold but, as the halo progenitors collapse, their asso- 
ciated baryons are shock-heated and begin to cool radiatively. In 
the bulge-dominated case, the baryon fraction decreases slightly at 
first but it soon settles into a quasi-steady state, with a cold frac- 
tion (stars plus dense gas with pg ^ 7 x lO^^'^gr cm~^) that is 
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Figure 3. Left: the evolution of the rms radius (top panel) and the specific angular momentum (middle and bottom panels) of the dark matter particles that lie 
within ry at 2 = 0. In the bottom panel the dashed line indicates L oc a?'!'^ and the dotted Une L = const. Right: as the left but for the particles that make the 
inner dark matter halo as defined in Section 2. All physical quantities are normalized to their values at the present day. 



approximately 2.5 times larger than the warm/hot fraction (the re- 
maining baryons). The total baryon fraction in the collapsed regions 
is in excellent agreement with the mean value of ~ 90% of the 
cosmic me an found in non-ra diative gas simulations of a ACDM 
universe bv lCrain et al. U2OO7I) . By contrast, in the disc-dominated 
case, strong winds eject almost half of the baryons from the pro- 
genitor halos by 2 ~ 4. Some of these, however, are recaptured 
after 2 ~ 1.5 and the baryon fraction recovers to about 75% of 
the cosmic mean by the end. The warm/hot phase behaves simi- 
larly to the bulge-dominated case, but the cold component is much 
lower, particularly during the crucial period between 2 = 4 and 1 
when much of the merger activity is taking place. The cold baryon 
fraction increases slightly towards the end, reaching ~ 0.5 of the 
universal value at the end. In what follows, we shall show how these 
differences in the history of the baryon content influence the final 
morphology of the galaxies in a profound way. 



3 RESULTS 

We now calculate the evolution of the physical (i.e. not the comov- 
ing) angular momentum of the particles that make up each of the 
three components defined in Section 2 at the present day: the total 
halo, the inner halo and the galaxy. We identify the same particles at 
each output time in the simulation, find the position and velocity of 
their centre of mass and compute the physical angular momentum 
of each system, Lt, as 

Lt = ^ li = ^ mi{ri X Ui), (1) 

i i 

where rt and Vi denote, respectively, the position and velocity of 
particle i relative to the centre of mass of the system. The total 
specific angular momentum, L, is given by L = Lt/Mt, where 
Mt is the total mass of the system. 



3.1 Dark matter component 

First, we focus in the dark matter component. The evolution of the 
baryons has only a minor effect on the evolution of the halo and 
so the results are very similar for the bulge- and disc-dominated 
galaxies as already seen in Fig. (2] We show results only for the 
latter case. 

The top-left panel of Fig.[3]shows the evolution of the physical 
rms radius of the dark matter halo as a function of the scale factor, 
normalised to the present day rms radius. The evolution closely 
follows the spherical collapse model: the system expands, reaches 
a maximum radius at 1 + 2 ~ 3 and begins to collapse at 1 + 2 ~ 2 
to form a virialized system of size about half the radius at maxi- 
mum expansion. The middle-left panel shows the evolution of the 
magnitude of the specific angular momentum, L — of all the 
dark matter particles that lie within r-v at 2 = as a function of a, 
once again normalised to the present day value. During the rapidly 
expanding phase, (1 + 2 > 4), the system gains angular momentum 
in proportion to a'^^^, as first calculated from tidal torque theory by 
[white (198%. During the collapsing phase, (1 + 2 < 2), the mag- 
nitude of the angular momentum remains nearly constant. To show 
the dependence on the scale factor more clearly, in the bottom- 
left panel of the figure we scale the specific angular momentum by 
a'^^^. This plot clearly shows the two distinct evolutionary phases, 
indicated by the dashed (L oc a'^^^) and dotted (L = const) lines. 
Our resul ts agree well with previous analyzes of N-body simula- 
tions (e.g. IWhitel 19841 : ICatelan & TheunjI 19961) . 

The right panels of Figure [3] show the same quantities as the 
left panels but now for the particles that make up the inner dark 
matter halo (i.e. the 10% most bound mass) at 2 = 0. Since this 
subsystem has a higher overdensity than the halo as a whole, it 
reaches maximum expansion earlier and begins to collapse at 1 + 
2 ~ 5. The inner halo shrinks more rapidly than the halo as a 
whole and ends up with a radius that is only about 1 /5 of the size 
at maximum expansion. At 1 + 2 '--^ 5, the particles that will end up 
in the inner halo are not contained in a single object, but in several 
subclumps which subsequently merge. This is illustrated in Fig. |4] 
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which shows the spatial distribution of the inner halo particles at 
1 + 2: ~ 4 and 3.3. Many of the fragments merge between these 
two epochs and it is these mergers that determine the evolution of 
the rms size of the system during this period. 

The middle and bottom panels of Fig. [3] show that the evo- 
lution of the specific angular momentum of the inner halo is very 
different from that of the halo as a whole. Initially, during the ex- 
pansion phase, the system gains angular momentum in much the 
same way as the halo as a whole, roughly following the L oc a'^^'^ 
scaling. However, after maximum expansion, (1 -f 2; < 4), the in- 
ner halo rapidly looses most of its angular momentum. This be- 
haviour is due to the intense merging activity illustrated in Fig.|4] 
As fragments are drawn towards the centre by dynamical friction, 
the asymmetric distribution of dark matter produced by gravita- 
tional tides exerts a torque which transfers the orbital angular mo- 
mentum of the fragments to the outer halo. Each merger event is 
accompanied by a decline in the angular momentum and in this 
way 90% of the angular momentum is drained from the inner halo. 
During this phase, the angular momentum of this system declines 
roughly as a~^. This beha viour was fi r st not ed in the early cold 
dark matter simulations of iFrenk et al.l jl985l) . The redistribution 
of ang ular momentum during viria lization has been studied in de- 
tail bv lD'Onghia & Navarrol bOOTl) . 



3.2 Bulge-dominated galaxy 

We now perform the same analysis on the baryonic component of 
the bulge-dominated galaxy illustrated in the left panels of Fig. [T] 
The total stellar mass of this system is 6.75 x 10^^ Mq and the 
total gas mass is 7.5 X lO^h-^MQ. 

The evolution of the rms radius of this galaxy, plotted in the 
upper-left panel of Fig.|5] is remarkably similar to that of the inner 
dark matter halo. The system expands until 1 -f z ~ 5 and then 
collapses. The evolution of the magnitude of the specific angular 
momentum (lower-left panel) is also very similar to that of the in- 
ner dark matter halo. During the expanding phase, the system gains 
angular momentum in proportion to a'^^^, as expected from linear 
theory. After 1 + 2 — 4, the system rapidly loses its angular mo- 
mentum, just as the inner dark matter halo did, transferring it to the 
outer halo. The decline in the specific angular momentum of the 
galaxy during the collapsing phase is less abrupt than for the inner 
halo, scaling roughly a L oc a~^. 

The similarity in the behaviour of the bulge-dominated galaxy 
and the inner halo arises because even before the collapse phase, 
the cold baryonic material is distributed in a similar way to the 
mass of the inner halo, as may be seen in Fig. [2] Radiative cooling 
of gas is very efficient in dense subclumps at early times and, in 
this simulation, much of the gas cools and turns into stars inside 
the fragments that will later merge to make the inner halo. Very 
little gas is left over for later accretion and, as shown in Fig.|2l the 
cold gas fraction hardly changes between z — 4 and the present. 
Just as in the dark matter only case, the fragments, now containing 
a mixture of dark matter and cold baryons (stars and dense gas), 
merge together transferring their orbital angular momentum to the 
outer halo. With most of the final cold baryons already in place 
within the fragments at early times, the end result of the collapse 
and merger phase is a bulge-dominated object. 



3.3 Disc-dominated galaxy 

We now analyse the disc-dominated galaxy illustrated in the right- 
hand panels of Fig. [T] This galaxy has 4.17 x W^'' Mq in stars 



and 1.53 x lO^'^h'^ Mq in cold dense gas. The results are shown 
in the right hand panels of Fig. [5] 

The top-right panel of Fig. [5] shows that the baryonic mate- 
rial that will end up in the galaxy expands until 1 + 2 ~ 4 and 
then begins to collapse. This happens slightly later than for the 
bulge-dominated galaxy. The lower-right panel shows that, as for 
all other objects and components, the evolution of the angular mo- 
mentum during the expanding phase follows linear theory closely, 
with L cx a^^^. However, unlike the bulge-dominated galaxy, the 
angular momentum of the system remains nearly constant during 
the collapsing phase (1 + 2: < 3), indicating that the accreting gas 
conserves its angular momentum. This behaviour is reminiscent of 
that of the halo as a whole, illustrated in Fig.jS] Indeed, the lower- 
right panel of Fig.|5]is very similar to the bottom-left panel of Fig[3] 
The correspondence between the evolution of the angular mo- 
mentum of the baryonic material of the disc-dominated galaxy and 
of the material that makes up the halo follows from the similar- 
ity of their spatial distributions at the time of maximum expansion. 
As discussed in Section 2, mergers of subclumps in this simulation 
induce large-scale shocks in the gas which result in an increased 
star formation efficiency. This, in turn, produces strong feedback 
which is further promoted by the assumed top-heavy starburst IMF, 
keeping the gas out of the merging fragments. Since the merger 
rate is higher at high redshift, this strong feedback suppresses the 
early collapse of gas into small proto-galaxies. The process is ap- 
parent in Fig. [2I the wind-driven ejection of gas at 2 4 dras- 
tically reduces the cold gas fraction in progenitor halos which re- 
mains much smaller than in the bulge dominated galaxy. Thus, most 
of the baryons that will make up the final galaxy become decou- 
pled from the sub-clumps early on and accumulate in a hot gas 
reservoir in the main halo. This hot gas has similar specific angu- 
lar momentum than the halo and, as it looses its pressure support 
through radiative cooling, it accretes onto the galaxy conserving 
its angular momentum. The key to the formation of disc galaxies is 
hence the suppression of the early collapse of baryons into merging 
fragm e nts, as anticipated by m any authors 1 Navarro & Steinmetg 
200d: lOkamoto et al.l |2005|; ISommer-Larsen & Dolgov .200 H : 



Thacker & CouchmanlboOlHWeil, Eke. & Efstathioulll998l) . 

Finally, in Fig. [6] we show the evolution of the unnormalised 
specific angular momentum of the whole halo (solid line) and of 
the baryonic components of both galaxies, the bulge-dominated 
(dashed line) and the disc-dominated (dotted line) objects. This 
figure indicates that it is not only the shape of the time depen- 
dence of the angular momentum that is similar between the halo 
and the baryons of the disc-dominated galaxy, but, remarkably, also 
the value of the specific angular momentum of the two compo- 
nents. Thus, our simulation provides strong support for the clas- 
sical theory of disc formation whereby tidally torqued gas is ac- 
creted into the centre of the halo conserving its angular momentum 
CFall & Efs tathiou 1980; Mo et al. 1998). This is the assumption 
made in semi-analytic models of galaxy formation (e.g. lCole et al. I 
I2OOOI) . By contrast, the bulge-dominated galaxy ends up with only 
about 25% of the specific angular momentum of its halo. 



4 DISCUSSION AND CONCLUSIONS 

The inability to obtain realistic discs by the present day in sim- 
ulations from cold dark matter initial conditions has been as- 
cribed to various possible ca u ses: (i) a lack of numerical reso- 
lution jKaufmann et alj 1 20071 : iGovemato et al. Il200 4l2007l) , (ii) 
an artificial transfer of angular momentum from a cold disc to 
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Figure 4. Projected positions of the dark matter particles which will eventually end up in the inner halo at 2 = 0. The left panel shows the distribution at 
1 + 2; ~ 4 and the right panel at 1 + 2 ~ 3.3. 




Figure 5. Evolution of the rms radius (upper panel) and specific angular momentum (lower panel) with scale factor for the baryonic component. The left 
panels correspond to the bulge-dominated galaxy and the right panels to the disc-dominated galaxy. All physical quantities are normalized to their value at the 
present day. 



a hot surrounding atmosphere ( IOkamotoet"ai]|2003l) . (iii) insta- 
bilities in a heavy disc which, by causing frag ments to form, 
end u p transporting angular momentum inwards (Robertson et al.l 
|2004) . and (iv) the tra nsfer of angular from merging substruc- 
tures to the outer halo jNavarro & Bend 1 199 ll ; In avarro & White! 
ll994lNavarro. Frenk. & Whitelll995h . The first thr ee of these four 
consid erations are unimportant in our simulations. Okamoto et al.' 
( I2OO3I) carried out a convergence study which suggests that, so long 
as the cold and hot gas phases are dynamically decoupled, resolu- 
tion effects are not important in calculations the size of those we 
analyse here. Such decoupling also suppresses process (ii). The 
adoption of an effective multiphase equation-of-state stabilizes a 
disc against gravitational instability thus alleviating (iii). Process 
(iv) is inevitable and is the dominant form of angular momentum 
transport (at least for the dark matter) in our simulations. 



The series of N-body/SPH si mulations of gal a xy fo rmation 
in a CDM universe carried out by lOkamoto et alj bOOSi) led to 
the conclusion that galaxies with very different morphologies can 
form in the same dark matter halo depending on the details of 
the assumed star formation and feedback prescriptions. In this pa- 
per, we have investigated the evolution of the specific angular mo- 
mentum as a possible explanation for why relatively small differ- 
ences in the assumed baryonic processes can result in such differ- 
ent outcomes. We have analyzed the bulge-dominated and the disc- 
dominated^ galaxies simulated (from identical initial conditions) by 
lOkamotoe t al. (2005). We trace back all the dark matter particles 
that end up within the virial radius of the halo (r„), as well as those 
that end up in the central parts of the halo by virtue of being in 
the top 10 percentile of the binding energy distribution. We also 
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Figure 6. Specific angular momentum of the dark matter component (solid 
line) and of the baryonic component of the bulge-dominated galaxy (dashed 
line) and of the disc-dominated galaxy (dotted line). 



components are very similar throughout the entire history of the 
galaxy. By contrast, the specific angular momentum of the bulge- 
dominated galaxy is only about 25% of the halo value. 

Our results confirm that the key to the formation of disc galax- 
ies in the cold dark matter cosmology is the suppression of the early 
collapse of baryons into small, dense halos. By contrasting the dif- 
ferent evolutionary histories of the angular momentum in bulge- 
and disc-dominated galaxies, our analysis explicitly reveals how 
regulating the rate at which gas is supplied determines the final 
morphology of the galaxy. A disc galaxy results when, as seen in 
Fig.|2l a large fraction of the baryon content of the halo is ejected 
by winds at early times. The ejected baryons (some of which are 
recaptured later) are thus decoupled from the subhalos that are des- 
tined to merge and join a hot gas reservoir which is tidally torqued 
in the same way as the host halo and thus acquires the same net 
specific angular momentum. Our results provide strong support for 
the classical theory of disc formation bv lFall & Efstathiou (1983) 
and IMo et al.l ( 1 19981) in which tidal torques impart the same spe- 
cific angular momentum to the dark matter and the gas and angular 
momentum is conserved during disc formation. 



trace back the baryonic particles that end up as a galaxy, i.e. as cold 
baryons (star particles and dense gas) within 10% of at z — 0. 

The angular momentum of the dark matter halo evolves al- 
most identically in the two simulations. The specific angular mo- 
mentum of the halo as a whole grows initially as predicted by tidal 
torqiie theory in the linear re gime, in proportion to a'^^^ jWhitd 
ll984ICatelan & Theunsll993) . This phase lasts until maximum ex- 
pansion is reached. Thereafter, the angular momentum of the dark 
matter remains constant, as the system collapses to a virialized con- 
figuration. By contrast, the evolution of the most bound 10% of the 
dark matter is strikingly different. In the linear phase this material 
behaves just as the halo as a whole but, after maximum expansion, 
its specific angular momentum rapidly declines to ~ 10% of its 
maximum value. The decline occurs in episodes associated with the 
mergers that build up the central halo. Most of the angular momen- 
tum of these fragments is invested in their orbits and is transferred 
to the outer halo by tidal forces as the fragments sink by dynamical 
friction. 

The baryonic component evolves very differently in the two 
simulations. In essence, in the disc-dominated galaxy, the specific 
angular momentum of the baryons tracks that of the halo as a 
whole, while in the bulge-dominated galaxy, it tracks the evolution 
of the central halo material. The reason for the difference can be 
traced back to the different strengths of feedback in the two galax- 
ies at early times. In the bulge-dominated galaxy, feedback is weak 
and most of the baryons rapidly cool and condense into stars within 
sub-galactic fragments. As these halo fragments merge and give up 
their angular momentum to the outer halo, so do the baryon clumps 
within them. The result is a slowly-rotating central bulge. A small 
amount of residual gas rains in later on forming a small disc. In 
the disc-dominated galaxy, by contrast, feedback is strong at early 
times, the gas is reheated before it can make a substantial mass of 
stars and, instead, accumulates in an extended hot reservoir which 
acquires a specific angular momentum similar to that of the dark 
halo. By the time this gas is able to cool, much of the merger activ- 
ity has subsided and so the gas is able to dissipate and collapse into 
an centrifugally supported disc configuration conserving its angular 
momentum. Not only does the angular momentum of the baryonic 
material track the evolution of the angular momentum of the halo, 
but the actual values of the specific angular momentum of both 
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